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Most of the theoretical results on the kinematic amplification of small-scale magnetic fluctuations 
by turbulence have been confined to the model of white-noise-like (5-correlated in time) advecting 
turbulent velocity field. In this work, the statistics of the passive magnetic field in the diffusion-free 
regime are considered for the case when the advecting flow is finite-time correlated. A new method 
is developed that allows one to systematically construct the correlation-time expansion for statistical 
characteristics of the field such as its probability density function or the complete set of its moments. 
The expansion is valid provided the velocity correlation time is smaller than the characteristic growth 
time of the magnetic fluctuations. This expansion is carried out up to first order in the general case 
of a d-dimensional arbitrarily compressible advecting flow. The growth rates for all moments of 
the magnetic-field strength are derived. The effect of the first-order corrections due to the finite 
correlation time is to reduce these growth rates. It is shown that introducing a finite correlation 
time leads to the loss of the small-scale statistical universality, which was present in the limit of 
the (5-correlated velocity field. Namely, the shape of the velocity time-correlation profile and the 
large-scale spatial structure of the flow become important. The latter is a new effect, that implies, in 
particular, that the approximation of a locally-linear shear flow does not fully capture the effect of 
nonvanishing correlation time. Physical applications of this theory include the small-scale kinematic 
dynamo in the interstellar medium and protogalactic plasmas. 



I. INTRODUCTION 



The study of the statistics of magnetic fluctuations excited by a random Gaussian white-noise- like advecting velocity 
field, was pioneered by Kazantsev K|, and, in more recent times, has generated a considerable amount of research 
(see, e. g., Ref. [|| and references therein, as well as Refs. p|-|Io|]). While much attention has concentrated on resistive 
dynamo problems, most often for very large magnetic Prandtl numbers, it is well known that the fundamental 
Zeldovich's, or "stretch-twist-fold," mechanism of the magnetic-energy amplification (the so-called "fast dynamo") 
is active regardless of the presence of the resistive (diffusive) regularization [^l) . If the initial seed magnetic field is 
concentrated on the scales of the same order as the characteristic scales of the advecting velocity, the stretching and 
folding of the magnetic-field lines by the random flow leads to an exponential growth of the magnetic fluctuations at 
scales that decrease exponentially fast, until the diffusive scales are reached [[l2|]o^,|po|| . This scenario is common in 
astrophysical applications such as the turbulence in the interstellar medium or in the protogalaxy where the Prandtl 
number ranges from 10 14 to 10 22 , giving rise to 7 to 11 decades of small (subviscous) scales available to the magnetic 
fluctuations |j|,|l3|,[l0]] . In fact, the initial diffusion- free regime may well be the only one practically important in such 
applications as far as the kinematic approximation is concerned, since the nonlinear saturation effects are likely to set 
in before the diffusion scales are reached On the fundamental physical level, the diffusion-free regime, in which 
the magnetic-field lines are fully frozen into the flow, exhibits most clearly the underlying symmetry properties of the 
passive advection |sj . 

With a few notable exceptions (such as Refs. [|l5 16 |J7|), the dominant approach in the existing literature on the 
turbulent kinematic dynamo problem has been to study the statistics of passive magnetic fields advected by a flow 
u(t, x) whose two-time correlation function is approximated by a <5 function, (u(i)u(t')) oc S(t — t'). This white- 
noise property of the velocity greatly simplifies matters: the evolution equations for such statistical quantities as the 
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correlation functions and probability density function of the magnetic field can be derived in closed form and yield 
themselves to exact solution. 

In this paper, we relax the white-noise assumption and explore the effects that arise when a finite-time correlated 
velocity field is introduced. This immediately raises the level of difficulty associated with solving the statistical 
problem. Within the theoretical framework adopted here, the difficulty can be described in the following terms. In 
the zero-correlation-time approximation, one essentially has to deal with only one closed differential equation that 
fully determines the desired statistics. Allowing for a finite velocity correlation time leads to an infinite number 
of interlinked integro-differential equations involving time-history integrals. These equations form an infinite open 
hierarc hy th at formally constitutes the exact description of the problem (these matters are explained in more detail 



in Sec. II A). Solving this hierarchy in its entirety without additional assumptions appears to be an impossible task. 



The most obvious way to make progress is clearly to try a perturbative approach, i.e., to consider the kinematic 
dynamo problem with an advecting field whose correlation time is short, but finite. If the correlation time r c is 
assumed to be small, one can expect to be able to construct an expansion in the powers of r c (in what follows, we will 
frequently refer to it as the t expansion) and calculate corrections to the growth rates of the moments of the magnetic 
field. This is the program that we undertake here. 

We consider the one-point statistics of the passive magnetic field in the diffusion-free regime. In this context, 
the infinite hierarchy we have mentioned above interrelates the one-point probability density function (PDF) of the 
magnetic field and an infinite set of response functionals. These are averaged multiple functional derivatives of the 
magnetic field with respect to the velocity field and its gradients. We develop a functional expansion method that 
allows us to calculate successive terms in the r expansion and derive in a closed form a Fokkcr-Planck equation for 
the one-point PDF of the magnetic field. We limit ourselves to advancing the expansion one order beyond the zero- 
correlation-time approximation. The result is a set of corrections to the growth rates of all moments of the magnetic 
field. These corrections are negative, so the growth rates are reduced. 

The expansion is carried out assuming that the velocity correlation time is small and keeping the time integral 
of the velocity correlation function fixed. The latter constraint ensures that the dynamo growth rate remains finite 
when the correlation time vanishes. An alternative way, which is sometimes deemed preferable on physical grounds 
(see, e.g, Rcf. JT^| ), is to fix the total energy of the velocity field. Since a 5-correlated velocity field must necessarily 
possess infinite energy, fixing the energy at a finite value leads to vanishing of the growth rates when t c = 0. The 
relative ordering of the terms in the expansion is, however, the same, regardless of what is kept fixed, so the technical 
side of the expansion method is unaffected. 

Our expansion technique will be given detailed treatment in the body of this paper. Here, let us rather discuss the 
finite-correlation-time effects that can be distilled on the basis of our approach. As it turns out, a number of new 
interesting phenomena manifest themselves already at the level of the short-but-finite-correlation-time approximation. 

In the case of the ^-correlated advecting flow, the one-point statistics of the passive magnetic field are universal 
in the sense that they only depend on one small-scale property of the velocity: the time integral of the one-point 
correlation tensor of its gradients, / dt (Vu(i)Vu(O)) . The essential novelty in the case of finite correlation time is 
that this small-scale universality is lost on two accounts. 

First, the r expansion exhibits a sensitive dependence on the specific shape of the time-correlation profile of the 
velocity field (in recent literature, this was first explicitly pointed out by Boldyrev ijTsjf ; see also Refs. [^9|,^0)). Namely, 
multiple time integrals of products of velocity correlation functions enter the expressions for the expansion coefficients. 
Choosing different correlation profiles leads to order-one changes in the values of these coefficients. The root of this 
nonuniversality lies in the topolog y of t he vertex-correction diagrams that contribute to the orders higher than the 
zeroth in the r expansion (see Sec. HE ). 

Second, the first-order terms of the r expansion feature a part that arises from the fourth-order derivatives of 
the velocity correlation function, i.e., from the second derivatives of the velocity field. In the one-point statistical 
approach, this is the first manifestation of the more general tendency that introducing finite correlation times brings 
into play the large-scale structure of the velocity field. A related effect is the loss of Galilean invariance due to the 
fact that the expansion terms also depend on the actual energy of the velocity field, i.e., on the rms value of the 
sweeping velocity. Indeed, now that the trajectories of the fluid elements have a "memory" of themselves, which 
extends approximately one t c back in time, we should naturally expect that there will appear an effective "correlation 
length" of the velocity (in what regards the one-point statistics of the fields it advects) approximately equal to ut c . 
Therefore, the one-point statistics of the passive fields now depend not only on the instantaneous velocity difference 
between two fluid particles that meet at a given time (i.e., the velocity gradient at a point), but also on the velocity 
that swept them into place and on the variation of the velocity gradient over the correlation length. This appearance 
of first-order corrections due to the second derivatives of the flow is a new effect, which indicates, in particular, that 
the customary approximation used in the Batchelor regime, where the advecting velocity is assumed to be locally 
linear 21 22] ], is only justified for the (5-correlated-in-time advecting fields. 

Such arc the main qualitative consequences of introducing a finite-time-correlated velocity field into the kinematic 
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dynamo problem (or, in general, any passive-advection model). A few words are in order as to the quantitative impact 
of a finite correlation time on the dynamo action. As we have already mentioned, the effect of the first-order corrections 
is to reduce the growth rates of all moments of the magnetic field. Besides the nonuniversal dependence on the spatial 
and temporal structure of the velocity correlation function, the reduction depends in a universally calculable way 
on the usual set of parameters: the order of the moment, the dimension of space, and the degree of compressibility 
of the flow. The overall magnitude of this reductive effect is measured by the expansion parameter, whic h is o f 
the order of r c 7, where 7 is the growth rate of the magnetic energy. It is not hard to demonstrate (see Sec. [II D) 
that r c 7<i ~ (T"c/Teddy) 2 , where r e ddy is the "eddy-turnover" time of the advecting turbulent velocity field and d is 
the dimension of space. In a standard Kolmogorov-type turbulence setting, one would, of course, expect any such 
approximation to be valid at best marginally, since r c ~ r e ddy Astrophysical plasmas offer more variety in this respect, 
as their driving forces (typically supernova explosions) can, in fact, decorrelate faster than the turbulent eddies turn 
over [^J. In any event, the small- t c expansion does not offer much more than qualitative, or, at best, semiquantitative, 
information about the way the dynamo action is modified by the finiteness of the correlation time. It is, of course, 
clear that introducing a finite correlation time cannot altogether suppress the fast-dynamo mechanism Jl(||7j . On the 
other hand, our conclusion that some reduction of the growth rate should be expected, is corroborated by numerical 
evidence |2^J^^4| that suggests a reduction of about 40% to 50%. In fact, in Sec. [V, we offer a semiquantitative 
evaluation of the finite- t c correction to the growth rate which yields a reduction of approximately 40% in the three- 
dimensional case and for r c ~ r c ddy- Of course, this is at best just an indication of the well-behaved character of our 
expansion, rather than a truly solid quantitative confirmation of it. 

The literature on the r expansion and finite-correlation-time effects is not extensive. Kliatskin and Tatarskii |25| ] 
were the first to propose the hierarchy of equations for the response functionals as a starting point for a method of suc- 
cessive approximations as applied to the description of waves propagating in a medium with random inhomogeneities. 
Vainshtein |26| ] applied this method to the mean-field kinematic dynamo theory. The Kliatskin-Tatarskii method and 
its relation to our functional expansion method are discussed at the end of Sec. p Oj . Van Kampen and Terwiel j27p 
developed the so-called cumulant expansion method; van Kampen's review article also contains a good critical 
survey of other r-expansion schemes predating his work. His method was later applied in the kinematic-dynamo 
context by Knobloch p^j and Chandran pj| . Their treatment was Lagrangian and did not include any effects due to 
the explicit spatial dependence in the induction equation. Consequently, the nonuniversality of the r expansion with 
respect to the s patia l structure of the velocity correlator was not captured. The van Kampen method is discussed 
Parallel to our development of the functional expansion method, Boldyrev |l£ 



in detail in Sec. [ID 



proposed a r- 

expansion method that was based on the exact solution of the induction equation in the Lagrangian frame and offered 
a way to calculate the second moment of the magnetic field that elicited the nonuniversal character of the r expansion 
with respect to both temporal and spatial properties of the velocity correlation tensor. Molchanov, Ruzmaikin, and 
Sokoloff []l6| considered the statistics of the kinematic dynamo in a renovating flow using the formalism of infinite 
products of random matrices. (See also Ref. |l5| ] for the treatment of the kinematic mean-field dynamo in a renovating 
flow.) A version of their approach was later advanced by Gruzinov, Cowley, and Sudan Q. Considerable progress was 
achieved in a nonperturbative way by Chertkov et al. |29| , who studied the passive-scalar problem in two dimensions 
for arbitrary velocity correlation times. However, their method only works in the two-dimensional case. 

Thus, while we now seem to have a fairly good understanding of the structure of the r expansion and such qualitative 
features as the loss of the small-scale universality, an adequate nonperturbative theory of the kinematic dynamo and 
passive advection in finite-time-correlated turbulent velocity fields remains an open problem. 

This paper is organized in the following way. In Sec. [n| our functional expansion method is systematically developed 
on the example of the simplest available passive-advection problem: that of the Lagrangian passive vector in an 
incompressible flow. In this model, no explicit spatial dependence is present. In Sec. [IA, Sec. [IB, and Sec. IIC, we 



present a functional formalism that allows one to systematically construct successive terms in the r-expanded Fokker- 
Planck equation. The dependence of the expansion coefficients on the specific functio nal f orm of the velocity time 
correlation profile emerges. The expansion is carried out up to the first order in r c . In Sec. IID, our method is compared 
with the van Kampen cumulant expansion metho d [^9| . We ascertain that results obtained via the van Kampen 
method are consistent with our s. Finally, in Sec. HE, we discuss the underlying structure of the r expansion in 
diagrammatic terms. In Sec. HI, the general arbitrarily compressible space-dependent dynamo problem is solved with 
the aid of the functional expansion. At this level, the nonuniversality with respect to th e spatial structure of the 
velocity correlations, as well as the loss of Galilean invariance, become evident. In Sec. [II A, we explain the emergence 
of an infinite hierarchy of equations for the characteristic function and various averaged response functionals of the 
magnetic field in the passive dynamo problem with finite-time-correlated advecting flow. The hierarchy is advanced 
up to the emergence of the second-order response functions. In Sec. IIIB, we construct the r expansion up to first 
order in the correlation time, which leads to a closed equation for the characteristic function of the magnetic field. 
In Sec. Ill C , wc derive the Fokkcr-Planck equation for the one-point P DF of the magnetic-field strength valid to first 
order in the correlation time. The distribution is lognormal. In Sec. [IID, we calculate the rates of growth of all 
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moments of the magnetic field with (negative) first-order corrections. Finally, in Sec. IV, we give a semiquantitative 
argument that relates the expansion parameter to the ratio of the correlation and eddy-turnover times of the velocity 
field. We also evaluate the finite- r c reduction of the magnetic-energy growth rate in a model incompressible turbulence 
consisting of eddies of a fixed size. In Appendix |a|, we provide the basic relations that allow one to express the results 
we have obtained in the configuration space in terms of the spectral characteristics of the velocity field. Some of the 
more cumbersome technical details of the r expansion are exiled to Appendix [b|. 

II. THE GAUSSIAN FUNCTIONAL EXPANSION FORMALISM 

In this Section, we explain the Gaussian functional method for constructing the short-correlation-time expansion 
for passive advection problems. Working out such expansions for specific problems often involves a fair amount of 
algebra, which tends to obscure the otherwise transparent ideas behind them. In an attempt at the maximum possible 
clarity of exposition, we first consider a model that, while preserving most of the essential features of the passive- 
ad vection problems, offers much greater technical simplicity. Namely, let us consider the following stochastic equation 
in d dimensions: 

d t B l = a\B k . (1) 

All the fields involved explicitly depend on time only. The specific initial distribution of B l is not important for the 
derivation or the validity of the results below. Spatial isotropy is always assumed. The matrix field <xl(i) is Gaussian 
with zero mean and a given two-point correlation tensor: 

(4(t)ai(t'))=TH K (t-t'), (2) 
1$ = 6^6 a + a(5i6j + 5\5{), a = -l/(d + 1). 

These equations can be interpreted to describe the evolution of a passive magnetic field in a Lagrangian frame, where 
the Lagrangian advecting velocity field is Gaussian and incompressible, and the tensor a], is its gradient matrix. In 
the more general context of the theory of passive advection, equations ([[]) and (0) model the stochastic dynamics of 
a vector connecting two Lagrangian tracer particles in an ideal fluid. 

We assume that the temporal correlation function n(t — t') of cr|(t) has a certain characteristic width r c , i.e., the 
field 0^(4) possesses a correlation time r c . Our task in this section is to construct an expansion of the statistics of B 4 (t) 
in powers of r c , which is assumed to be small. The limit r c — > ought to be taken in such a way that the time integral 
of the correlation function is kept constant: 

°° K 

dr k(t) — — — const and t c k <C 1. (3) 

The white-noise limit of zero correlation time is realized by setting k(t) — kS(t). 

The first step in our averaging scheme is to define the characteristic function of the field B l (t) , 

Z(t\n) = (Z(t- lf i)) = (exp^B^)]). (4) 

Here and in what follows, the overtildes designate unaveraged random functions. Upon differentiating Z(t\ \x) with 
respect to time and making use of Eq. (|l|), we obtain a new stochastic equation: 

d t Z = aiix~Z = k\uiZ, (5) 

where the auxiliary operator A k has been introduced for the sake of notational compactness. 

Our objective now is to learn how to obtain a closed equation for the averaged characteristic function Z(t; fx), i.e., 
how to average Eq. (||) when a\ has a nonzero correlation time. The inverse Fourier transform (with respect to //,) 
of the resulting equation will be the Fokker-Planck equation for the PDF Pit; B) of the passive field B % . 

A. The Hierarchy of Response Functions 

We start the construction of the functional expansion by developing an exact formalism that describes the one-point 
statistics of the field B l (t). Let us average both sides of Eq. (||) and "split" the mixed average that arises on the 
right-hand side with the aid of the well-known Furutsu-Novikov (or "Gaussian- integration" ) formula |30|,[ni : 
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d t Z(t) = Al(ai(t)Z(t)) = k\Tf a \ At' «(t - f)G%(t\t% (6) 

Jo 

where we have suppressed the /i's in the arguments, used the formula for the second-order correlation tensor 
of cr|(i), and introduced the averaged first-order response function: 

W-<a?MO>-(^). m 

This function is subject to the causality constraint: Gg(t\t') = if t' > t [hence the upper integration limit in Eq. (||)]. 
Integrating Eq. (||) from to t, taking the functional derivative 5/5o~P(t') of both sides, averaging, setting t' = t, and 
taking causality into account, we get 

G%(t\t) = A%Z(t). (8) 

We have thus obtained the equal-time form of Gp(t\t'). In order to determine the response function at t 1 < t, we 
simply take the functional derivative 6/6o~P(t') of both sides of Eq. (0) and find that each element of the unaveraged 
tensor G^(t\t') satisfies an equation identical in form to Eq. (^). Upon averaging this, we obtain an evolution equation 
for Gp(t\t') subject to the initial condition (||) at t — t': 

d t G1(t\t')=A l j (ai(t)G^t\t')). (9) 

This equation can now be handled in the same fashion as Eq. (||), the average on the right-hand side split via 
the Furutsu-Novikov formula in terms of the correlation tensor of of (t) and the appropriately defined second-order 
averaged response function G**? 2 (t\t\, *2). At equal times, the latter can be expressed in terms of G%(t\t') just 
as Gp(t\t) was expressed in terms of Z(t). At different times, we obtain the evolution equation for GV^(t\ti,t2) by 

taking the functional derivative of the equation for Gg(t\t') and averaging. 

An infinite linked hierarchy can be constructed by further iterating this procedure and introducing response functions 
of ascending orders. Let us give the general form of this hierarchy. Define the nth-order averaged response function: 

G ^"" ^( fcftfrffUfr.) )- (10) 

This function has two essential properties: (i) it is causal: G^ 1 "'p n (t\ti . . . t n ) = if any U > t; (ii) it remains invariant 
under all simultaneous permutations of the times t±, . . . , t n and indices a±, . . . , a n , @±, . . . , (3 ni which correspond to 
changes of the order of functional differentiation in the definition (|l(]) . The nth-order response function satisfies the 
following recursive relations: if t\, . . . , t n < t, 

dtG a p[::Z •••*«) = k i T it + + \ f d ^+i ^ - Wi)^ 1 m • ■ ■ Wi); (n) 

J o 

if, say, t n = t and ti, . . . , t n -i < t n , 

Gpl'.'.'.Pn-lPn ^l* 1 • ■ ' = k p"G fll'.'.'.Pnll ' ' ' ( 12 ) 

The hierarchy is "forward" at different times and "backward" at equal times. The characteristic function Z(t) is 
formally treated as the zeroth-order response function. 



B. The White-Noise Approximation 

The white-noise approximation is obtained by setting n{t — t') = iiS(t — t'). We are then left with just Eq. (^), where 
the time history integral reduces to ^K,Gg(t\t), which is substituted from Eq. (||). This produces a closed evolution 
equation for Z(t). The Fourier transform of it is the Fokker-Planck equation for the PDF of B l at time t in the 
(^-correlated regime: 

dtP(t) = \ Ttk\k%P{t) = \ LP(t). (13) 
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In order to be not overly burdened by notation, we typically use the same symbol for denoting an operator in the 
Fourier space of the /it's and its analog in the configuration space of the B's. This should lead to no confusion, as the 
context will always be clear. Thus, 



dB 



■B k = -{ 5. 



dB 1 



Due to isotropy, the probability density function P depends on the absolute value B 
in Eq. (fL3|) can therefore be written in the following isotropic form: 



f _ jii/3 i k i a _ 



d-l 



B 2 



dB 2 



(14) 

|B| only. The operator L 
(15) 



which turns Eq. ( |l3| ) into the familiar Fokker-Planck equation for the one-point PDF of the magnetic field in the 
kinematic ^-correlated dynamo problem taken for the incompressible flow [^|j9| . The resulting distribution is lognorma! 
and the moments of B satisfy 



dt (B ri 



I d-l 



2 d+ 1 



n (n + d) R (B n ) 



(16) 



the Lagrangian and Eulerian statistics are the same 



This is the expected outcome, because, as was shown in Ref. 
for the (5-correlated incompressible flow. 

Thus, the solution in the (5-correlated limit is quite elementary. Things become much more complicated once the 
white-noise assumption is relaxed and a nonzero, however small, velocity correlation time is introduced. 



C. The Recursive Expansion 

In order to construct an expansion in small correlation time, it is convenient to combine the equations ( |TT| ) and (JlJ 
into one recursive integral relation that expresses the nth-order response function in terms of its immediate precursor 
and its immediate successor: 

+ ^ T tt + + \ f d*' T dWi «(f - tn^Gf^llWW ■ ..*»+!)• (17) 

Jtn JO 

The above relation is exact and valid for t\,... ,i„-i < t n < t. Due to the permutation symmetry of the response 
functions, this does not limit the generality. The desired expansion is constructed by repeated application of the 
formula (17). 

Let us substitute the formula ( |l7j ) with n = 1 for the first-order response function into the right-hand side of Eq. (|TJ) : 



d t Z{t) = L dti K(t - ti)Z(ti) 
Jo 

+ T^TZ^^kn fdh fdt' f dt 2 K(t-h)K(t' -t 2 )G%%(t%,t 2 ), (18) 
JO Jt\ Jo 

where the operator L is defined in (|l5|). We now use the formula ( |l7| ) to express the second-order response function 
on the right-hand side of the above equation: for t 2 > t\, we have 

G?^'\hM) = ^lG%{h\t l ) 

+ llT^l f dt" f dh K (t" - t 3 )G a p ^(t"\t u t 2 ,t 3 ), (19) 
Jt 2 Jo 

while for t 2 < t\ we flip the variables, t\ ^> t 2 , to make sure that the first-order response function on the right-hand 
side do not vanish: 



G 



Gft^(f|*i,*») = A^(tiN 

+ A|If£ f At" f dt 3 n(t" - ts)G%%% (t"\t 2 , t lt t3). (20) 
Jti Jo 

The recursion relation (|l7|) is now applied to the first-order response functions in the formulas ([l9|) and ( pp| ) : 

G^(t2|ti) = A«^(* 1 )+A«T^| f^dt" f dt 3K (t"-t 3 )G^(t"\ tl ,t 3 ), (21) 

Jti Jo 

G^{t x \t 2 ) = kf % Z{h)+klT^l r&t" f dt 3K (t"~t 3 )G a ^(t"\t 2l t 3 ). (22) 

Jt 2 Jo 

All this must be substituted into Eq. (p"§|): 

d t Z(i) = L / dii n(t - ti)Z(ii) 
Jo 

+ ^ 1 1 ^ 2A " A ™ A ft A ft Ati /*df f 1 dt 2K (t- tl )K(t'-t 2 )Z(t 2 ) 

Jo Jti Jo 

+ r&T^%2&k a £k% fdh fdt f dhKit-uW-hWh) 







J2(t), (23) 



where the remainder i?(£) contains the assembled terms that involve quintuple time integrals. 

So far, all the manipulations we have carried out have been exact. It is now not hard to perceive the emerging 
contours of the small- r c expansion. Since the time-correlation function n(t — t') is a profile of width ~ r c , the area 
under which is constant and equal to R, the triple time integrals in Eq. (|2^) are of the order of t c R 2 , while the quintuple 
time integrals absorbed into R(t) are of the order of t 2 R 3 . Further application of the recursion formula (JTt]) to the 



second- and third-order response functions in the equations (|19|)-(22) leads to the appearance of more multiple time 
integrals of the time-correlation function n(t~t'). These integrals are of orders t 2 R 3 , t 3 R 4 , etc. We would only like to 
keep terms up to first order in the correlation time. The remainder term R(t) in Eq. (E3) can therefore be dropped. 

Remark on the physics of the t expansion. The above argument is based on the stipulation made at the beginning 
of this Section that the r expansion must be carried out keeping the integral of the velocity time correlation function 
constant [formula (g)]. This r equir ement is natural because it leads to finite dynamo growth rates in the limit of 
zero correlation time [see Sec. |n§ Eq. (H|)]. However, it is also acceptable to institute an alternative, arguably 



more physical, requirement that the total energy of the velocity field (i.e., the rms velocity) remain constant (as, 
e.g, in the numerics of Ref. fl7j| ). Quantitatively, this means that k(0), rather than J ° o dTK(r), is kept fixed. Un- 
der this constraint, the terms that we have previously estimated to be of orders R, t c R 2 , t 2 R ; , etc., and hence, R 
being constant, to represent the zeroth, first, second, etc. orders of the r expansion, should now be reevaluated as 
follows. Since R ~ t c k(0), these terms are of orders t c k(0), t 3 k(0) 2 , t^k(Q) 3 , etc., and therefore constitute the first, 
third, fifth, etc. orders of the expansion. The shortcoming of this approach is that the dynamo growth rates vanish 
when r c = 0, so formally there is no nontrivial zero-correlation-time limit. With R — const, this problem was avoided 
because the energy was formally infinite when r c = (a 5-correlated velocity field cannot have a finite energy). In any 
event, we see that, since the difference between keeping R and k(0) constant does not affect the relative magnitudes 
of the terms in the expansion, our expansion scheme remains valid in both cases. Let us therefore proceed with our 
construction. 



The dependence of the right-hand side of Eq. (|23| ) on the "past" values of Z (i.e., on its values at times preceding t) 



can also be resolved in the framework of the small- r c expansion. Formally integrating Eq. (23), we get, at times t\ < t, 



Z(h) = Z(t)- L f dt' f dt 2 K(t' - h)Z(h) + ■ ■ ■ (24) 



ti Jo 



Upon substituting this onto the right-hand side of Eq. (|23|) and again discarding all the terms of orders higher than 
the first in r c , we get 
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d t Z{i) = L dtiK(t-h)Z(t) 
Jo 

ft rt ft' 



dt x \ dt' dt 2 K(t-t 1 )K(t' -t 2 )Z(i) 



'0 Jt x JO 

pt pt pt\ 

+ (L 2 - Li) dh dt' dt 2 n(t-t 1 )K{t' -t 2 )Z(t) 

JO Jt x Jo 

+ (L 2 -L 2 )f dh f dt' f dt 2 K(t-ti)K(t' -t 2 )Z{t). 



Jo Jt 

We have introduced the following two operators: 



L2 = ^ i TZ\[kf i ,L]=2di 



a 2 

02 



d+l 



(25) 

(26) 
(27) 



where the square brackets denote commutators. We see that the terms in Eq. ( p5| ) that contain L 2 cancel out, and 
only those terms remain that are due to the non-self-commuting nature of the operator A*. 

Finally, we inverse- Fourier transform Eq. ( p5|) into the B space and take the long-time limit, t ^> t c . The following 
Fokker-Planck equation with constant coefficients results: 



LP, 



where the coefficients. 



K 2 



— ^ lim / dh / dt 2 [ dt 3 K(t — ti)K(t 2 — t 3 ), 
r c K 2 *->°° Jo Jt! Jo 

ft ct P t 2 



— rrr lim 

T C K Z t— >oo 



dh / dt 



dt 3 K(t - tt)K(t 2 - t 3 ), 



(28) 

(29) 
(30) 



are constants that depend on the particular shape of the time-correlation function n{t — t'). Thus, the r expansion is 
nonuniversal in the sense that the specific choice of the functional form of the small-time regularization directly affects 
the values of the expansion coefficients (cf. Ref. (l^]). As an example, let us give the values of the coefficients K\ 
and K 2 for two popular choices of n(t — t'): 



K (t - *') = 5— exp[-|i - t'\/r c ] => K 1= K 2 = 0.5 



K(t - t') 



exp[-(f - t') 2 /T 2 } Ki w 0.33, K 2 w 0.23. 



(31) 
(32) 



In Sec. Ill, we will apply the method we have presented above to the more realistic general compressible kinematic 



dynamo problem in the Eulerian frame. 

Remark on the Kliatskin-Tatarskii method. The Gaussian hierarchy given by the equations ( |TT| ) and ([l2]) and 
based on repeated application of the Furutsu-Novikov formula was proposed by Kliatskin and Tatarskii |^5j as a basis 
for constructing successive-approximation solutions of the problem of light propagation in a medium with randomly 
distributed inhomogeneities. Their method in its original form was carried over to the mean-field dynamo theory 
with finite-time-correlated velocity field by Vainshtein [p6| . The method we have outlined in this section, while 
also based on the response- function hierarchy (|Tl|)-(|l2|), differs substantially from that developed and applied by 
these authors. Their successive-approximation scheme consisted essentially in writing out the first n equations in the 
hierarchy ([ll])-(|l2]) and then truncating it at the nth step by replacing n(t — t n +i) by RS(t — i„+i) in the equation for 
the nth-order response function. This gave a closed system of equations that could be solved. Carried out in the first 
order, such a procedure would correspond to setting n(t' — t 2 ) — RS(t' — t 2 ) in Eq. ( |l8|) and consequently n(t 2 — t 3 ) = 
R8(t 2 — t 3 ) in the expressions ( |29] ) and ( |30| ) for the coefficients K\ and K 2 . Such a substitution leads to K\ = 
and K 2 = (2/t c k) J q dr tk(t), which is incorrect. The reason for this discrepancy is that, in the time integrals 
involving multiple products of the correlation functions n(t — t\), n{t 2 — t 3 ), etc., the latter cannot be approximated 
by S functions plus first-order corrections even in the small- r c limit. 
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D. Comparison with the Van Kampen Cumulant Expansion Method 



The evolution equation (^|) for the "unaveraged characteristic function" Z(t; (a) is a stochastic linear differential 
equation whose form agrees exactly with that of the general such equation considered by van Kampen jl9 20 and 
simultaneously by Terwiel p^] : 



dtZ(t) = A(t)z(t), 



(33) 



where, in our case, A(t) = A^cr|(t). In his work, van Kampen developed a formalism that allowed one to construct 
successive terms in the short-correlation-time expansion of Z = (Z) in terms of the cumulants of the operator A. 
Terwiel's projection-operator method was shown by its author to be equivalent to that of van Kampen. 

Let us see what happens if the passive-advection problem given by Eq. (|l|), or, equivalently, by Eq. (||) is subjected 
to van Kampen's expansion algorithm. The latter proceeds as follows. 

Start by writing the formal solution of Eq. (0) in terms of the time-ordered exponential: 



Z(t) = 



exp / di'A(i') 



Z(0) 



1+ / dt 1 A(t 1 )+ / dii / dt 2 A(ti)A(t 2 ) + 



Z(0). 



(34) 



This solution is averaged assuming that the initial distribution of Z is independent of the statistics of A: 



Z(t) 



dti 



di- ; 



dii 

*3 



dis 



dt 2 (A(ii)A(i 2 )> 

dt4(l(*l)i(*3)A(i 3 )i(*4)) 



Z(0). 



(35) 



Here all the odd-order averages have vanished (recall that A = Afcr^,). The closed equation for Z(t) is now obtained 
as follows. First, the formal solution (pq) is differentiated with respect to time: 



dt 1 



d t Z(t) 



dU 



dh (A{t)A(h)) 
dt 3 {A{t)A{t{)A{t 2 )A(tz) 



Z(0). 



(36) 



Second, Z(0) is expressed in terms of Z(t) by formally inverting the operator series on the right-hand side of Eq. (|35|), 
whereupon Z(0) is substituted into Eq. (J36]). Keeping only the terms that contain up to three time integrations, as 
we did in the previous section, we get 



d t Z(t) 



uo 



d*i (A(t)A(ti)) 



dii 



dU 



dii / di 2 
'o 



di 3 (A^Ait^A^Aih)) 
di 3 {A(t)A{t 1 )){A{t 2 )A(u)) + 



Z(t). 



(37) 



The quadruple average in the above expression splits into three products of second-order averages in the usual Gaussian 
way. Since 



we have 



(A^Aih)) = n(t - ii) T^A^Aj = n(t - h) L, 



(A(t)A( tl )A(t 2 )A(t 3 )) = n(t - h)K(t 2 - t 3 )T%T£ n AfAjA^A' 
+ n(t i 2 ) K (ii - i 3 ) TfrTfr A^A^Al 



(38) 
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+ K(t 


- h)n{ti 


- h)Tt q T? p m A^A^Al 


= n{t 


- ti)n(t 2 


-t 3 )P 


+ K(t 


- h)n(ti 


- t 3 ) (L 2 - L x ) 


+ n(t 


- h)n{ti 


-t 2 )(L 2 -L 2 ), 


= n{t 


- h)n(t2 


-h)L\ 



(39) 
(40) 

The operators L, Li, and L 2 are the same as those in the previous section [see definitions (p"5|), (p6|), and (P7|)]. The 
averages (|38|), (39), and ( p0| ) are now substituted into Eq. (|37|). The triple time integrals can be argued to represent 
(all of the) first-order terms in the small- r c expansion in the same way as it was done in Sec. II C. In the limit t 3> t c , 
the coefficients in Eq. ([37]) do not depend on time. The resulting Fokker-Planck equation for the PDF P(t;B) [which 
is the inverse Fourier transform of Z(t; fi)] obtained by the van Kampen method and analogous to Eq. (E8h is then 



dtP 



1 



Ci + C 2 J - X - t c k(C -Ci- C 2 )L 



LP, 



2d+l 

where the coefficients Co, Ci, and C 2 , that depend on the shape function n(t — t'), are as follows: 





4 








/"*2 




C = 




lim 




/ dt 2 


/ dt 3 «(t - 


- tl)K(*2 




T C K 2 


t—>oc 




hi j 







4 












Ci = 




lim 




r dt 2 


/ d£ 3 k(* 


- t 2 )K(ti 




T C K 2 


t—>oo 




'o 


Jo 




4 












c 2 = 


-9 


lim 


(* 


r dt 2 


/ dt 3 k(* 


- h)n{tx 






t— -*oo 




h 


Jo 



(41) 

(42) 
(43) 
(44) 



By comparing the definition of Co with those of the coefficients K\ and K 2 in Sec. II C [see formulas ( |29| ) and (30)], we 
immediately establish that Co = K\ + K 2 . Furthermore, it is also not hard to ascertain that Ci = K\ and C 2 = K 2 . 
Therefore, the last term in Eq. ( |4l| ) vanishes, and the first-order Fokker-Planck equations ( ]2q) and (^) are identical. 
Thus, the results obtained via the van Kampen method are consistent with ours. Unlike the van Kampen method, 
however, our method does not involve any nontrivial operator algebra and is therefore better suited for a wide variety 
of applications. In particular, the stochastic equations containing spatial derivatives (such as the convective derivatives 
present in all Eulerian passive-advection problems) can be handled without much additional difficulty (this will be 
done in detail for the full kinematic dynamo problem in Sec. III). 



E. Discussion: The Vertex Corrections 



While the particular methods one employs to obtain the successive terms in the r expansion may vary and depend 
on one's taste and the specific demands of the stochastic problem at hand, the underlying structure of the r expansion 
remains the same and is rooted in the common properties of all turbulence closure problems (see, e.g., Ref. |3^| ). As we 
have stated in general ter ms in the intro duction to this paper, and as was clear from our const ructi on of the response- 
function formalism in Sec. II A -Sec. II C or of van Kampen's explicit series solution (35) in Sec. IIP , averaged solutions 
of stochastic equations such as Eq. (£3|) can be represented in terms of infinite sums of multiple time-history integrals 
containing products of time-correlation functions ft(tj — tj) in the integrands. This summation can be visualized in 
terms of Feynman-style diagrams. The n-point diagrams represent the terms containing n time-history integrations. 
As an example, Fig. [j] lists the three possible four-point diagrams. 

It was noted by Kazantsev [Q (see also Ref. [0) that the white- noise approximation corresponds to the partial 
summation of all ladder-type diagrams such as the four-point one shown in Fig. [l](a) . The distinctive property of these 
diagrams is that the pairs of points U, tj at which the time-correlation functions in the integrands of the time- history 
integrals are taken, can be fused without interfering with each other. No essential information is therefore lost when 
the time-correlation functions n(ti — tj) are approximated by S functions. However, in all orders of the r expansion 
but the zeroth, diagrams with more tangled topology appear: e.g., in the first order, these are the diagrams [j](b) 
and 0(c)]. Such diagrams are often referred to as the vertex corrections. Fusing points in these diagrams leads to the 
loss of terms that cannot be neglected [ p3[ . This is the context in which the emerging nonuniversality with respect to 
the shape of the time-correlation profile should be viewed. 

In this paper, we restrict our consideration to the first-order terms in the r expansion. The relevant diagrams are 
the four-point ones shown in Fig. |[ The diagrams [|(b) and [|(c) give rise to the coefficients Ci and C 2 , respectively 
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[see Eq. (41) and formulas 
£3 - * t± in the diagram 
and formulas (E9T) and 



43| ) and (E3)]. Upon changing variables t\ <-> £2 in the diagram [l](b) and — > £2, £2 — > fa 
cl , we see that these diagrams equally well correspond to the coefficients K\ and K% [Eq. (|2g 






FIG. 1. The three possible fourth-order diagrams. The times ti, £2, t$ are integration variables and may float along the axis, 
but their positions relative to each other and to the time t are determined by the limits of integration and therefore preserved. 
A dashed line connecting any two times ti and tj represents the time-correlation function fc(t, — tj) in the integrand of a triple 
time-history integral. 



III. THE FUNCTIONAL EXPANSION FOR THE KINEMATIC DYNAMO IN A 
FINITE-TIME-CORRELATED VELOCITY FIELD 

In this Section, we use the functional expansion method developed in Sec. || to construct the r expansion for the 
general diffusion-free kinematic dynamo problem in the Eulerian frame with an arbitrarily compressible velocity field. 

Through the convective derivative, an explicit spatial dependence is now present in the problem. This leads to 
the appearance of the new effect advertised in the Introduction: while the zeroth-order terms in the expansion only 
depend on the one-point correlation properties of the velocity gradients, the first-order terms also depend on the 
energy of the advecting velocity field and on the one-point correlation function of its second derivatives. The former 
represents the loss of Galilean invariance, the latter the loss of the small-scale universality and the advent of the 
sensitive dependence of the statistics on the large-scale structure of the velocity correlations. 

In this Section, all statistics are Eulerian. For the questions regarding the transformation of PDF's of passive 
fields from the Eulerian to the Lagrangian frame, we address the reader to Ref. ||, as well as to Ref. p8|, where the 
r expansion is treated as a problem in stochastic calculus and Lagrangian statistics are discussed. 



A. The Gaussian Hierarchy 

The magnetic field passively advected by the velocity field u l (t, x) evolves according to the Hertz induction equation 
(formally in d dimensions): 

d t B l = -u k B\ + u\B k - u k k B\ (45) 

where u\. = du l /dx k , B\ = dB l /dx k , and the Einstein summation convention is used throughout. Let the advecting 
velocity field u l (i, x) be a homogeneous and isotropic Gaussian random field whose statistics are defined by its second- 
order correlation tensor: 

xV(i',x')> = K ij (t-t', x-x'), (46) 
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where, as a function of the time separation t — t' , the correlator is assumed to have a finite width r c , which we 
will call the velocity correlation time. As we will only study the one-point statistics of the magnetic field, all relevant 
information about the velocity correlation properties is contained in the Taylor expansion of K y around the origin: 

« ij (r, y) = ko(t)*« - \ «2(r) [y 2 S^ + 2ayV] 

+ ^K 4 (T)y 2 [y 2 S^ +2byW] +•••, (47) 

as y — > 0. Here a and & are the compressibility parameters. Between the purely incompressible and the purely 
irrotational cases, they vary in the intervals 

We should like to mention here that the choice of the coefficients of the small-scale expansion (^) of the velocity 
correlation tensor is, strictly speaking, not entirely unconstrained. As k u (t, y) is a correlation function, it must 
be an inverse Fourier transform of a proper correlation function in the Fourier space |34| . In Appendix |A|, we give 
the expressions for the coefficients of the expansion (^) in terms of the spectral characteristics of the velocity field. 
We further note that, while a and b can certainly be functions of t, we will not overly shrink the limits of physical 
generality by assuming that they are either constant or slowly-varying functions of time, i.e. that they do not change 
appreciably over one correlation time. 

In order to determine the one-point statistics of the magnetic field, we follow the standard procedure fjEf and 
introduce the characteristic function of J3 l (i,x) at an arbitrary fixed spatial point x: 

Z{t;ji) = (i(t,x;/x)) = (exp^B^x)]). (49) 

As usual, the angle brackets denote ensemble averages and the overtildes mark unaveraged quantities. The function Z 
is the Fourier transform of the PDF of the vector elements B % . Due to spatial homogeneity, Z does not depend on 
the point x, where B i (t,x.) is taken. Upon differentiating Z with respect to time and using Eq. ([45|), we get 

d d 

d t Z = -u k Z M + u\ni—Z - u k kfM — Z = -u k Z M + A k u l k Z, (50) 



where, for the sake of future convenience, we introduce the operator 



d d 

Af = Mi^ S k fii — , (51) 



which will turn up repeatedly in this calculation. 

In order to obtain an evolution equation for the characteristic function Z(t; /u) of the random magnetic field, we 
average both sides of Eq. (|50|). Since, due to the homogeneity of the problem, {u k Z^) = —(u k k Z), we may write the 
equation for Z(t; fi) in the following form: 

d t Z(t)=(5? + A k )(u%(t,x)Z(t,x)) 

= -(S k +A k ) J\ tl Jd^^it-h^-^G^it,^,^), (52) 

where the mixed average on the right-hand side has been "split" with the aid of the Furutsu-Novikov ("Gaussian- 
integration" ) formula |3^,[nJ , and the fj, dependence in the arguments has been suppressed for the sake of notational 
compactness. We have introduced the first-order averaged response function of the following species: 



G£(i,x|*i,xi) = (G^ftxlti.xi)) = ( ff; X) ) . (53) 



As a response function, satisfies the causality constraint: (t, x|ti, xi) = for t\ > t. The same-time form 
of Go* can be obtained in terms of the characteristic function Z(t): integrating Eq. ( |5C| ) from to t\, taking the 
functional derivative 5/8u^ l (£', xi), averaging, setting t' = t±, and taking causality into account, we get 
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G^ i 1 (ti,x|ti,xi)=<y(x-x 1 )A° 1 Z(ti). (54) 

In order to find G^(t,x\t\,Xi) at t > ti, we take the functional derivative S/Su 13 ^ it' , x') of both sides of Eq. ( |50| ) 
and establish that each element of the unaveraged tensor Ga* satisfies an equation identical in form to Eq. (JHT 



dt&fc (t, x|i 1)Xl ) = - M m (i, x)G^ im (t, x\ti, x x ) + A>£(t, x)^ 1 (t, x|t 1)Xl ). (55) 

Subscripts such as " lTn " in the above equation mean, in accordance with the usual notation, the partial differentiation 
with respect to x m , viz., d/dx m . 

We must now average Eq. (|55| ) in its turn, to obtain an evolution equation for G^(t, x|ti,xi) at t > t\. Using the 
initial condition (]54|), let us write this evolution equation in the integral form valid for all t > t\\ 

G%(t,x\t 1 ,x 1 ) = S(x-x 1 )A%Z(t 1 ) 

- [ dt' [ dt 2 /d d x 2 L m& (i , -t 2 ,x-x 2 )G^ 1 /32m (t',x|i 1 ,x 1 ;t 2 ,x 2 ) 
Jti Jo J L 

+ (*' - * - x 2 )G^ 2 2 (f , x|t!, x i; i 2 , x 2 )] , (56) 

where the mixed averages have again been "split" by the Furutsu-Novikov formula, at the price of introducing two 
new second-order response functions: 

G p 1 p 2 ( t ,x|ti,x i; i 2 ,x 2 ) = I —j— — — > , (57) 

G2 1 ? a (t / ,x|t 1 ,x 1 ;i2,xa) = ( — 5 L^jl^I \ . ( 58 ) 

01 " 2 \^ 1 (t 1 ,x 1 )^ 2 fe,x 2 )/ 

In the same way that the equal-time first-order response function was expressed in terms of Z(t) [Eq. (ESjJ)], the 
second-order response functions at t' = t\ or t' = t 2 can be expressed in terms of G^ {t\, x|i 2 , x 2 ) or G^ 1 (t 2 , xjii, xi), 
respectively. Because of causality, the former representation would be valid provided t\ > t 2 , the latter in the 
opposite case t 2 >t\. At other times, ti,t 2 < t', the functions G^ 1 ^ and G^ 1 ^ 2 satisfy integral equations analogous 
to Eq. d56|), where third-order response functions make their appearance. An infinite open hierarchy can thus be 
obtained by further iterating this procedure and introducing response functions of ascending orders. This hierarchy 
constitutes the exact description of the statistics of the kinematic dynamo problem with arbitrary velocity correlation 
time. 



B. The t Expansion 

The expansion in small correlation time must be carried out in such a way that the time integral of the velocity 
correlator k^(t, y) remains constant. Since K L i(T,y) has a finite (small) width r c , we can conclude that the double 
time integral on the right-hand side of Eq. ( ]56] ) must be of first order in the correlation time r c . As we are only 
interested in constructing the r expansion up to first order, it is now sufficient to calculate the second-order response 
functions G^ 1 ^ and G^ 1 ^ 2 with zeroth-order precision. 



We have already mentioned that recursive relations completely analogous to the relation (56) can be derived for the 
second-order response functions. The latter are thereby expressed as their equal-time values plus double time integrals 
of the same sort as that which appeared on the right-hand side of Eq. (|56|). These time integrals are first order in 
the correlation time and can therefore be neglected. The equal-time values of the second-order response functions 
are obtained by formally integrating Eq. (pOj), taking functional derivatives of it, averaging, and using causality. The 
second-order response functions are thus expressed to zeroth order in terms of the first-order ones. These latter can 
by the same token be replaced by their equal-time values, which only contain the characteristic function Z(t). The 
resulting expressions, valid to zeroth order, must be substituted into the first-order term (the double time integral) 
in Eq. (|56|). All these manipulations, which require a fair amount of algebra, are relegated to Appendix^. Here we 
simply give the resulting expression for the first-order response function, valid to first order in r c : 
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G^(f,x|t 1 ,x 1 ) = <5(x-x 1 ) 











+/' 


,ti 

dt' / dt 2 K 




Jo 


<9 2 <5(x - xi) 


f dt' f dt 


Q x mQ x n 


Jt x Jtt 



'"(*' -t 2 ,0)A^Z(t 1 ) + O(r c 2 ). (59) 

This expression must now be substituted into the time-history integral on the right-hand side of Eq. (^|[). This 
gives a closed integro-differential equation for the characteristic function Z(t). However, the dependence on the past 
values of Z is spurious and can be resolved to first order in r c . Indeed, we can formally integrate Eq. (|52| ) from t\ 
to t and, using the zeroth-order value of the first-order response function [the first term in the formula (p9|)1, get 

Z(t 1 )=Z(t) + (6^ +A£) fdt' f dt 2 n^{t'^t 2 ,Q)k^Z{t 2 ) + 0{r 2 c ). (60) 

Jti Jo 

The double time integral in this equation is of first order in t c , as usual. 

Upon assembling the equations (|||), (^9]), and (|6C|), we finally arrive at the following closed partial differential 
equation for Z(t): 



Jo 

dii f dt' [ 'diaK^ (*-*!, 0) 
Jti Jo 

x (f - t 2 , 0) (5* + A* ) [A^ 1 , A» ] A£ Z(t) 

/ dh f dt' f dt 2 Kf ai (t-h,0) 
Jo Jti Jti 

x - t a ,0)(tf + A?) [A-, + A£)A£]Z(t) 



ti 



dt x / dt' I dt 2K f ai (t-t u 0) K %%(t'-t 2 ,0)(5* + A*)A%Z(t) 

Jtx Jo 

- / dh fdt' [ dt 2K % imn (t-t 1 ,0)^"(t'-t 2 ,0)(^ + A^)A^Z(t). (61) 



The square brackets denote commutators. 

Note that, besides the second derivatives of the velocity correlation tensor, the first-order terms contain the fourth 
ones, as well as the undifferentiated tensor itself [in the last term in Eq. fl6l|)]. The latter implies the loss of Galilean 
invariance, the former the loss of small-scale universality in the sense that the large-scale structure of the velocity 
correlator starts to play a role. This effect could not have been captured if the velocity field had been assumed to be 
purely a combination of the instantaneous velocity at a given point and a linear shear. 



C. The Fokker-Planck Equation 



In order to obtain the Fokker-Planck equation for the PDF of the magnetic field, we must inverse-Fourier trans- 
form Eq. ( fill ) back to B dependence. The inverse Fourier transform of Z(t; fi) is the one-point PDF P(t; B). We will 
continue using the symbol A^ to denote the counterpart of the operator A^ in the B space: 

A^d-ljtf-B^ + tfB'^. (62) 
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Due to the isotropy of the problem, the PDF P(t; B) will in fact be a scalar function of the field strength B only. 
Thus, all the operators that appear on the right-hand side of the B-space counterpart of Eq. ([n]) must, after they are 
convolved with the velocity correlation tensors, be expressible in terms of B. Let us use the Taylor expansion ( Ejj ) of 
the velocity correlator to calculate the tensor convolutions in Eq. (^l|). We have 



K ij (r,0) = ko(t) S ij , 
k%(t, 0) = -n 2 (r) [5^d kl + a (SiS* + 
; j Wmro (r,0) = M(r)[2(d + 



- K 2(r)T^, 
b )^S kl + (d + 4)b(6lSi+SlSi)] 



(63) 
(64) 

(65) 



A number of second-order differential operators (with respect to B) arise in Eq. (|6l]). In the zeroth-order term, we 
have 



£ = + A?) A* 




(l + P)B-^ + (d+l)(3 



(66) 



two operators appearing in the first-order terms result from the non-self-commuting nature of the operator A k [see 
the second and the third terms in Eq. (pj|)] p6|: 



rpi/3i rpmf3 2 ( Xk , \k\ [ \ai An]ia 2 
1 ka 1 1 na 2 \°i + J H ) V V f3i ' 1V m\ ly j3 2 



d 2 (d- 1) 
(d+l) 2 

rpifil rpmf3o 



2d(d- 1) 
d+l 



d 



B 



dB 



1 + ^J B ^ 



^ 2 = T^r„T 2 2 + A*) [A* , (C + A" ) A£] = T*£ + A?) [A£ , L] 



d 2 



(67) 



(68) 



and, finally, there are two other operators due to the presence of the convective term (i.e., explicit spatial dependence) 
in the induction equation [see the fourth and the fifth terms in Eq. (J6l|)] : 



d(d-l) 
(d + l) 2 



B 



d_ 

dB 



2/3 d(d + 1)-1 



1 + ^ 



•AfUS 1 



d 3 



a (d + i) 2 



2(d-l)(d + 4) / d 



01 

d 2 + 4d+S 
2d(d + 4) 



d + 3 



d 2 
d 



(69) 



C B 



dB 



B dB 

(d + 2)(d + 3) 
2(d + 4) 



(70) 



In all of the above, /3 = d[l + (d + l)a] and C, = d[2 + (d + 3)6] are compressibility parameters that vanish in the case 
of incompressible flow. In this latter case, the operators defined above simplify considerably: 



L\ = 

My = 



d+l 
d 



L, L 2 = 2dL, 

— L, M 2 = 2 ^ d + 4 h. 
1 ' d + 3 



(71) 
(72) 



If we take the long-time limit, i.e., t ^> r c , the coefficients in Eq. ( pi] ) do not depend on time t. We can now use 
the inverse Fourier transform of Eq. (|6|) taken in this limit and the isotropic operators listed above to assemble the 
Fokker-Planck equation for the PDF of the magnetic field. This equation contains the desired corrections that are of 
first order in the velocity correlation time r c and represent the first available manifestation of the finite-correlation-time 
effects. We have 



d t P 



11 
2 



L- 1 - t c k 2 [K x (Lx + Mi) + K 2 L 2 + K 2 M 2 ] ) P, 



(73) 
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where the overall dimensional factor is 



R 2 = 2 / drK 2 (r), 
Jo 



and the coefficients 



K 2 
K 2 



4 
4 



' c ^2 

4 



T C K2 t^co 



dt 3 n 2 (t 


- ti)n 2 (t 2 


-*3), 


dt 3 k 2 (t 


- ti)n 2 (t 2 


-ta), 


dt 3 Ki(t 


- ti)no(h 





(74) 

(75) 
(76) 
(77) 



are constants that depend on the particular shapes of the time-correlation functions kq(t), k 2 (t), and Ka[t). Such 
sensitive dependence is a new feature and represents a loss of universality with respect to the specific time-correlation 
profiles (cf. Ref. fL8|). As we have pointed out in Sec. [II B, the universality with respect to the functional form of 
the velocity correlator in space is also lost (this effect is incorporated into the coefficient K 2 ). 

Let us also list the much more compact form that the Fokker-Planck equation ( |73| ) assumes in the case of an 
incompressible velocity field: 



d t P 



11 
2 



<«>^ iA mr k > 



LP. 



(78) 



Here it is especially manifest that the true expansion parameter in the problem is T c R 2 d. This is a general statement 
that holds regardless of the degree of compressibility, as can be readily verified by counting powers of d in the general 
expressions for the operators L, L%, L 2 , Mi, and M 2 [formulas (|6q)-(|70|)1. 

It is evident that the distribution resulting from Eq. ( |73|) is lognormal, which is a well-known fact in the kinematic- 
dynamo and passive-advection theory. Since we are interested in the quantitative description of the fast-dynamo 
effect, we will now proceed to calculate the growth rates of the moments of the magnetic field. 



D. The Dynamo Growth Rates 

The evolution of all moments of B can be determined from Eq. (|75| ) . The nth moment is calculated according to 

{Bn) = T{d/2) I ABBU+ P ^B). (79) 
Upon multiplying both sides of Eq. ( |73| ) by B d+n ~ 1 and integrating over B, we find that (B n ) satisfies: 
9t (B n ) = 7 (n) (B n ) 

= y {r(n) - r c S 2 d[Jfiri(n) + K 2 T 2 (n) + K 2 T 2 (n)] } (B n ) , (80) 
where the nondimensionalized zeroth-order growth rates are (cf. Ref. B) 



T(n) = ^±n[n + d+(n-l)f3]., 



(81) 



and the universal parts of the (negative) first-order corrections arising from the second- and fourth-order terms in the 
velocity correlator (|47|) are 



ri(n) 
r 2 («) 

f 2 (n) 



1 d- 1 



2 d+ 1 
d-l , n ( 8 
—I n{n + d) { 1+ dl 



(d- l)(d + 4) 
d{d + 3) 



, n ( d 2 + Ad + 2 



(82) 
(83) 
(84) 
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We observe that, for n = 0, T = Ti — T2 = f 2 = 0. This simply means that both zeroth- and first-order terms in the 
t expansion preserve the normalization of the PDF, i.e., our expansion is conservative, as it should be. 

In the incompressible flow, the total growth rate of the nth moment can be written in a more compact form: 



K2 d — 1 . . 

^ n) = TdTT n{n + d) 



~ 1 2 1 2 d(d + 3) 



(85) 



We see that the corrections to the growth rates of the magnetic-field moments are negative, so the growth rates are 
reduced. The amount of reduction depends on a variety of factors including the dimension of space, the order of the 
moment, the degree of compressibility, the functional form of the velocity correlator in time and space, and, of course, 
the velocity correlation time. Let us note that our general results derived for an arbitrarily compressible velocity 
field reveal no qualitatively essential effect of compressibility on the behavior of the first-order finite-correlation-time 
corrections to the dynamo growth rates in the diffusion-free regime. Compressibility of the flow simply leads to 
additive (and positive) corrections to the incompressible values of T(n), Ti(n), r 2 (n), and r 2 (n). Quantitatively, 
these corrections may affect the exact conditions for the break-down of the first-order approximation. For more 
discussion of the compressibility effects in the kinematic dynamo (with a ^-correlated velocity field), we address the 
reader to Refs. PH||3"8[ . 

We remind the reader that here we have studied magnetic fluctuations in the diffusion-free regime and therefore 
dropped the term in the induction equation that is responsible for the resistive regularization. Such an approach is 
justified for plasmas with very large magnetic Prandtl numbers (e.g., the ISM or the prototogalaxy) and applies to 
the initial stage of the small-scale dynamo that lasts for a time of order t ~ log Pr that elapses before the magnetic 



fluctuations reach resistive scales 14 10. After that, or if the Prandtl number is of order unity or small (as is, e.g. 
the case for the Sun), resistive effects must be taken into account. In this case, the calculation of the moments of the 
magnetic field via the Fokkcr-Planck equation for its PDF as presented in this Section does not apply because of the 
closure problem associated with the diffusion term [the equations for Z(t,x;fi) and Z(t;n) do not close]. However, 
the general r-expansion method proposed in this paper can, in principle, be applied to multipoint correlators of the 
magnetic field, for which treating the diffusive case presents no conceptual difficulty. One-point moments can then 
be obtained by fusing the points at which the multipoint correlators are taken (cf. Refs. [ p9pC| ]). Although it is the 
diffusive case that is studied in most numerical simulations, where Pr rarely exceeds 100, it is not necessarily the 
most relevant one in the context of the (proto)galactic dynamo, for which Pr ~ 10 14 -r 10 22 . Indeed, as we already 
pointed out in the Introduction, the initial (proto)galactic seed field may well be strong enough for the kinematic 
approximation to break down while the dynamo is still in the diffusion-free stage If this is the case, the effect of 
magnetic diffusion must be studied in conjunction with nonlinear saturation of the magnetic fluctuations fll7|. 



IV. A PHYSICAL EXAMPLE: THE ONE-EDDY MODEL 



In real astrophysical environments, such as the interstellar medium and the protogalactic plasmas, the magnetic 
fields are acted upon by a Kolmogorov-like turbulence with a fully developed inertial range about three decades 
wide (Re ~ 10 4 ). While the velocities of the turbulent eddies excited by the Kolmogorov cascade decrease with the 
scale of the eddy, the velocity gradients increase (see, e.g., Ref. pl|). Therefore, the dominant role in the process 
of amplification of the small-scale magnetic fluctuations is played by the smallest eddies. With this circumstance 
in mind, one often considers, for modeling purposes, a synthetic incompressible turbulent velocity field consisting of 
eddies all of which have the same fixed size but random isotropic orientation (for detailed discussions of the galactic 
and protogalactic dynamo, we refer the reader to Refs. ||[l3 14 lC|j) . In this Section, we will present a brief discussion 



of the implications of the r-expansion theory developed in Sec. for such a model problem, which will henceforth 
be referred to as the one-eddy model. 

The velocity field in the one-eddy model is specified as follows 

d d k 



u%x) = J -^e ik -V(i,k), (86) 
where the Fourier modes u l (t, k) are random variables that satisfy 

(u\t, k)u j (t', k')) = (27r) d ,5(k + k') (s ij - 5(k - k ) K (t - t'). (87) 

In this case, k%(t) cx k(t), and, upon using the relations listed in Appendix |a|, we get 
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1 (d-l)(d+2) 
Ko(T) = fcf d+1 K2(T) ' 



, 2 d + 3 



K4(T)=fc ° 2(d + 4)(d+l) K2(t) ' (89) 



Let us specify a plausible velocity time-correlation profile: 

K2 (r) = || exp (-M) . (90) 

For this correlation function, which corresponds, for example, to the well-known Ornstein-Uhlenbeck random process 
(see, e.g., Ref. 0), the coefficients of the r expansion ( plf ) are ifi = 7^2 = 1/2. The relations ( |S8] ) and ( p9[ ) provide 
the value of K 2 . 



(d-l)(d + 2)(d+3) 
2(d+l) 2 (d + 4) 



^ = ^StP ^- (91) 



Let us define the "eddy-turnover" time r e ddy ~ (^o u ) 1 of such a velocity field according to the following relation: 
' ..iir^-Ol-^-^™^ (92) 



'eddy Tc J -™ T c 1 



where we have used Eq. (|88|) to express ko in terms of k 2 . Note that the same expression is obtained if r c ddy ~ (Vu : 
V u )-!/ 2 i s formally defined in terms of the velocity gradients (without recourse to the one-eddy model): 

1 - 1 ^drl^-OI- jg^M * (93, 



T 2 7 , 

'eddy ' c ./-co 



We recall that the zeroth-order growth rate 70 of the magnetic-fluctuation energy (-B 2 ) is [see formula (p5f) 



Formulas (B3) and (p4) then imply 



7"cddy 



(d-l)(d + 2)_ tntS 
d+1 ^ (94) 



(d-l)(d+2) _ 

= T c 7 d = —— T c n 2 d. (95) 

d+1 



We have established a correspondence between the small parameter that has arisen in our expansion of the dynamo 
growth rates and the "physical" small parameter, which is the ratio of the correlation and eddy-turnover times. Of 
course, the above expression hinges on the definitions (j92] ) or ( J93] ) of r e ddy A simple physical argument can be made 
in favor of these definitions and the resulting formula (Pq). Namely, let us observe that when r c ~ r e ddy the eddy only 
stretches the magnetic field line in one of the d available directions during one turnover time, whence t c j ~ 1/d. The 
same estimate follows from the formula (j95[). 

Let us now evaluate the first-order correction to the growth rate of the magnetic energy. In the one-eddy model, 
one gets, upon using formulas ( |85| ) and j9l|) and taking K\ = K 2 = 1/2 for the Ornstein-Uhlenbeck time-correlation 
profile @, 

7 = 7W = 70(1-^), Ci = J*±J^. (96) 

We note that in three dimensions, Cd — 23/60 ~ 40%. When r c ~ r ddy, we have r c jod ~ 1, and the resulting 
growth-rate reduction of ~ 40% is in a good qualitative agreement with the available numerical results |^3|,[l7],|4|] . Of 
course, as we have already stressed in the Introduction, our r expansion is not designed for the case of r c ~ r c dd y , so 
the fact that it gives a fairly reasonable prediction should not be considered as an adequate quantitative corroboration 
of our theory. At best, one might conclude that the first-order expansion is well behaved for not-too-small values of 
the expansion parameter. 

Let us emphasize, however, that such a well-behaved expression has resulted from a number of essentially arbitrary 
(albeit physically reasonable) specifications of the parameters involved in the r expansion. One of the most physically 
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important points that we have tried to make in this work is, in fact, that the inclusion of finite-correlation-time 
effects leads to nonuniversal statistics, so the quantitative predictions of the theory can and will change appreciably 
if such factors as the shapes of the time-correlation profiles are changed. Namely, one would obtain expressions 
of the form (|9q ) with different values of the coefficient Cd- For sufficiently large values of T c 7od, the validity of the 
expansion (Bq) will break down, and the expression in the brackets may even become negative. However, the following 
heuristic argument can be envisioned in this context. 

Let us recall that the finite-correlation-timc effect was due to the presence of time-history integrals such as those 
that appear in the equations ( |52| ) , (|5^) , and (|5|) . The first-order corrections in the Fokker-Planck equation (|73| ) arose 
from systematically approximating the time evolution of the statistical quantities [response functions and characteristic 
function Z(t;/jL)] that entered these time-history integrals. The corrected ("true") value of 7 represents, in a rough 
way, the rate at which these quantities change. It would appear then that a better estimate of 7 would be obtained 
if 70 in the first-order term in the brackets in Eq. ( |96| ) were replaced with the corrected value 7. With this caveat, 
we would find that 



' l + C d r clo d' y ' 

To first order, this formula is equally accurate as Eq. (|96|). However, it better represents the fact that, as T c jod in- 
creases, the corrected value of of the growth rate should be expected to saturate [jl7|]. Of course, such considerations 
cannot substitute for an adequate nonperturbative theory of the passive advection and kinematic dynamo in finitc- 
time-correlated flows, which remains an open problem. 
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APPENDIX A: SMALL-SCALE-EXPANSION COEFFICIENTS OF THE VELOCITY CORRELATION 

TENSOR IN TERMS OF VELOCITY SPECTRA 

In this Appendix, we list the basic formulas that relate the coefficients of the small-scale expansion © of the 
velocity correlation tensor to the spectral characteristics of the velocity field. These relations allow one to apply the 
results on the small- r c expansion obtained in Sec. Ill to velocity fields that are specified in the Fourier, rather than 
configuration, space. They also provide a set of consistency constraints that must be respected when the specific 
functional forms of kq(t), «2(t), and Ki(r) are chosen. 

Let the advecting velocity field be given as a sum of spatial Fourier modes, 



d d k 



and let the Fourier coefficients w l (t,k) be random variables that satisfy 



{u l (t,k)u:>(t'X)) = {2n) d 6(k + k') 



K{k,t-t')S ij + H{k,t-t') 



k 2 



(Al) 



(A2) 



For the incompressible flows, k(k,r) = — k(/c,t); for the irrotational ones, K(k,r) = 0. The coefficients of the 
expansion (J47|) can then be expressed as follows: 



k 2 (t) 



d d k 
W) 



-g [dn{k,T) +h(k,T)], 



d d k 



d(d + 2) J (2ir) d 



k 2 [(d+2)K{k,T)+k(k,T)]. 



(A3) 
(A4) 
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1 



d d k 



2d(d + 2){d + 4) J {2ir) d LV J y ' ' y ' lj ' 



o(t) = k 2 (t) 1 



d d fc 



d(rf + 2) J (2ir) d 



k 2 k(k, t), 



d d fc 



6(r) = K4(rr d(d+ 2 )(d+4) y (2^f fc T) ' 

where the k-space integrations of radial functions can, of course, be written more explicitly as 



d d k 



(2ir) d (2ir) d J Q 



dkk 



d-l 



s d = 



27T rf / 2 

r(d/2)' 



(A5) 
(A6) 
(A7) 

(A8) 



The derivation of the above relations is straightforward and based on the expressions for the correlation functions 
of isotropic fields in configuration space in terms o f th eir spectra. For the 3-D case, these expressions can be found 
in Ref. A detailed derivation of the formulas ( A3)-([A7|) for the d-dimensional case is also given in Appendix A 
of Ref. toff. 



APPENDIX B: SECOND-ORDER RESPONSE FUNCTIONS 



In t his Ap pendix, we provide the zeroth-order expressions for the second-order response functions that we used 
in Sec. 1IIB. They are all derived in the same fashion: Eq. (pfl) is formally integrated, functional derivatives of it are 
taken with respect to the velocity field u l or its gradients u l k at the appropriate moments, the result is averaged, and 
the causality property of the response functions is used. Here we simply list the results. 

When t 2 > t%, we have, to zeroth order in the correlation time t c , 



G tip2 (*', x|*i , xi ; t 2 , x 2 ) = G^ 2 (t 2 , x\t i , xi ; t 2 , x 2 ) 



•5(x-x 2 ) 



A^G^xl^xO, 



= -5(x-x 2 )G^ /32 (t 2 ,x|t 1 ,x 1 ) 

Gf«l (t\ x|t x , xi; t 2 , x 2 ) = Gf«l (t 2 , x| t! , Xl ; t 2 , x 2 ) 
= -A" 2 (x - x 2 )G^ >/32 (t 2j x|ti , xi) + *(x - x 2 )A« 2 G^ (i 2) x|ti, Xl ), 

where we have introduced the following notation: by definition, 

5u m (t,x) 



Su a 2 (t 2 ,X 2 ) 



S%S(t-t 2 )A a *( x - X2 ) 



(Bl) 
(B2) 

(B3) 



The function A Q2 (x — x 2 ) is nonrandom and has the following property, which will be all that we need to know 
about it: 



^A« 2 (x-x 2 ) = <S«M(x-x 2 ). 



(B4) 



When t\ > t 2 , the expressions (Bl) and (B2) vanish by causality, so we have to flip the order of functional 
differentiation: 



G 7ifo (*', x|*i, xi ; t 2) x 2 ) = Gp"^ (t i , x|i 2 , x 2 ; 1 1 , x x ) 
= -A" 1 (x - xi)G fe!/3l (ti, x|i 2 , x a ) + S(x - xjfcfcGfr (t x , x|i 2 , x 2 ), 

G££ (t*, x|t! , xi ; t a , x 2 ) = G££ fa , x|t 2 , x 2 ; t x , x x ) 
= -A" 1 (x - xi)G« 2 !ft (t x , x|f a , x 2 ) + 5(x - x x )A^G^ (*i, x|t 2 , x 2 ). 



In Eq. (B5), the following obvious notation was used: 

G /32ft(* 1,X l <2 ' X2;il ' Xl ) = 



5 2 Z(ti,x) 



fe /32 (i 2 ,x 2 )^ Q 1 1 (<i,x 1 ) 



(B5) 
(B6) 

(B7) 



20 



and a new first-order response function appeared: 

G/? 2 (*i>x|i 2 ,x 2 ) 
The equal-time form of this function is 

G / 3 2 (t 2 ,x|<2,X 2 ) 



Su^(t 2 ,X. 2 ) 



dx 1 



■£(x-x 2 ) 



(B8) 



(B9) 



The first-order response functions that appear in the formulas (Bl), (|B2|), (|B5|), and (B6) can be written as their 
equal-time values (54) and (By) plus first-order terms. To zeroth order, we have therefore: for t 2 > t\, 



G Tifc (*'> x l*i) x i;*2,x 2 ) 



G^ 2 2 (i / ,x|i 1 ,x 1 ;t 2 ,x 2 ) 



d 



dx^ 



<5(x - xi] 



5(x- xi) 
.9 



<5(x - xi; 



5(x-x 2 )A° l Z(ti) 
A Q2 (x-x 2 )A^ 1 Z(i 1 ) 



5(x-x 1 )5(x-x 2 )A^A« i 1 Z(i 1 ); 



(BIO) 



(Bll) 



for tj, > t 2 , 



G^ 1 1 /Sa (t / ,x|ti,xi;t a ,x 2 ) 



-A Q1 (x-xx) 



6(x - x 2 ) 



A^(* 2 ) 



5(x-xi^ 



<9x" 



<5(x - x 2 ) 



A£A^(f 2 ), 



G^ 2 2 (i', x|t!, Xi ; * 2 , x 2 ) = -A" 1 (x - Xl ) 



5 



0aA 



<5(x- x 2 ) 



A^(* 2 ) 



+ *(x-xi)*(x-x 2 )A^A^Z(t 2 ). 



(B12) 



(B13) 



These expressions must be substituted into Eq. (|5fj). The volume integrals with respect to x 2 can be done, taking 
into account the extremely useful fact that all odd spatial derivatives of the velocity correlator k*- 7 (t, y) vanish at the 
origin (at y = 0). The results are: for t 2 > ti, 



J d d x 2 (t' - t a , x - x a )G£ ft m (f, x|t!, x i; f 2 , x 2 ) 



9 2 (5(x- Xl ) 



■ K 



mf3 2 (+1 



(f - 1 2 , o) + <$(x - Xl ) <;f„ 2 (<' - 1 2) o) A 



A^(*i), 



for ti>t 2 , 



J d d x 2 k™% (i' - f 2 , x - x a )G££ (^xl*!, x i; t 2> x 2 ) 
= <5(x - Xl ) <*(t' - t 2s 0) A^A^Zfa); 

| d rf x 2 K mfe (t' - 1 2 , x - x 2 )G^ 2)m (f', x|t!, x i; f 2l x 2 ) 
<5(x - Xl ) (f - i 2 , 0) + - *». 0) A^ 1 ] A^(i 2 ), 

/ d rf x 2 (f' - ta,x - x 2 )G^ 2 2 (i', x|t!, x i; fa, x 2 ) 
= <5(x - x x ) (f - t 2 ,0) A™ A^A^Z(ia). 



(B14) 



(B15) 



(B16) 



(B17) 



With the aid of these expressions and Eq. (pq), one obtains Eq. (p9|) of Sec. [II B 
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